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Summary 


The  frequency  response  of  a  shell  and  tube  heat  ev  ih^nger  was 
obtained  by  the  solution  of  an  ordinary  differential  equ'-'f ion  in  the 
frequency  domain.  This  equation  was  obtained  by  a  Four  Or  transformation 
of  the  set  of  partial  differential  equations  represent' ig  the  dynamic 
time  response  of  the  system. 

Three  simple  transfer  functions  were  evaluated  ;or  their  ability  to 
approximate  the  dynamics  of  the  exchanger.  The  trnrv'jr  function 

(1  -  Ke~tOS) 

(TjS  +  1)  (t2S  +  1) 

was  found  to  approximate  the  system  dynamics  exceptionally  well. 

A  method  of  evaluating  the  parameters  of  the  approximate  transfer 
function  by  relating  the  frequency  response  to  the  time  response  via 
Parraval's  theorem  was  also  evaluated.  Little  advantage  was  four.a  in  this 


INTRODUCTION 


Often  in  process  simulation  studies  and  process  control  it  is 
desirable  to  have  an  approximate  transfer  function  that  could  be  used 
for  dynamic  simulation  via  analog  computer  or  for  feedforward 
control  algorithms. 

In  many  cases  it  is  only  desired  to  know  the  system  response  to 
some  disturbance  at  a  fixed  space  corrdinate,  even  though  the  system 
is  space  dependent.  For  example,  if  the  system  under  consideration  is 
a  plug  flow  reactor,  the  reaction  occurs  continuously  along  the  length 
of  the  reactor  causing  reactant  and  product  concentrations  to  be  space 
dependent.  Perhaps  it  is  necessary  to  know  how  changes  in  the  feed 
composition  affect  product  distribution.  This  being  the  case  an 
approximate  transfer  function  relating  product  distribution  dynamics 
to  changes  in  feed  composition  could  be  employed  in  a  feedforward 
control  strategy. 

Sya terns  whose  dynamic  response  are  space  dependent  are  fee- 
quently  referred  to  as  distributive  parameter  systems.  The  equations 
describing  their  dynamic  behavior  are  partial  differential  equations 
in  time  and  one  or  more  space  coordinates.  Classical  time  domain 
solutions  for  these  equations  are  very  hard  to  come  by,  and  ofte'; 
they  are  considered  a  hopeless  undertaking. 

The  scope  of  this  work  is  concerned  with  developing  a  general 
analytical  technique  that  can  be  easily  implemented  in  obtaining 
approximate  transfer  functions  for  distributive  parameter  systems  of 


the  convective  type,  where  the  spacial  dependence  of  the  dependent 
variable  is  of  little  interest  from  a  control  viewpoint. 

This  technique  relies  on  a  distributed  parameter  model  Chat  is 
transformed  into  ordinary  differential  equations  in  the  frequency  domain. 
These  equations  are  then  analyzed  in  the  frequency  domain  using  well 
known  classical  techniques.  The  frequency  solution  of  the  ordinary 
differential  equation  provides  the  necessary  data  to  carry  out 
regressions  on  the  parameters  of  the  postulated  transfer  function. 

After  an  acceptable  fit  is  obtained  in  one  or  more  of  the  transfer 
functions,  they  are  inverted  into  the  time  domain  and  compared 
against  the  numerical  solution  of  the  model  equation. 
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Basic  Theory 


The  equations  describing  the  dynamics  of  distributive  parameter 
systems  are  partial  differential  equations  in  time  and  one  or  more 
space  coordinates.  Classical  time  domain  solutions  of  these  equations 
have  been  found  only  for  a  few  cases  of  practical  interest. 

It  has  been  known  for  many  years  that  linear  partial  differential 
equations  can  be  analyzed  in  the  frequency  domain  provided  the  system 
dynamics  are  described  at  a  fixed  position  cf  the  spacial  coordinates. 
Although  restricting  the  analysis  to  a  fix  point  in  the  spacial 
coordinates  presents  no  serious  drawbacks,  the  complex  frequency  ex¬ 
pressions  that  result  in  most  applications  are  not  suitable  for  time 
domain  inversion  with  classical  methods.  The  ever  Increasing  availa¬ 
bility  of  digital  computers  coupled  to  the  development  of  sophisticated 
non-linear  regressing  techniques  have  provided  the  necessary  tools  to 
make  the  time  domain  inversion  possible. 

In  this  section  the  theoretical  aspecte  involved  in  the.  fre¬ 
quency  analysis  of  the  partial  differential  equation  and  the  time 
domain  inversion  of  the  resulting  frequency  expression  will  be  discussed. 

To  translate  the  partial  differential  equation  to. the  frequency 
domain,  it  is  necessary  to  obrain  the  Laplace  transform  of  the  indi¬ 
vidual  terms  in  the  equation. 

The  Laplace  transform  of  a  continuous  and  differentiable  space 
dependent  variable  is  defined  as: 

CO 
n 

**  j  £(k,c)  exp  (wts)  dfc  (1? 

o 


F(x,s) 


By  using  equation  (1)  to  obtain  the  Laplace  transform  of 

5f one  obtains ; 

3t 

.  T  exp(.ts)dt  (2) 

Integrating  equation  (2)  by  parts  gives: 

CO  1  OD  CO 

Bf  (x,  t)  exp(_tg),jt  =  f  (x,t)exp(-ts)  I  +  s  i  f  (x,  t)exp(-st)dt  (3) 
J  Ot  'O  i! 

o  o 

Evaluating  limits  and  using  equation  (1)  to  obtain: 


=  sF(x>s)-f(x,o)  (4) 

The  Laplace  transform  of  higher  order  partial  derivatives  of 

fch 

f(x,t)  with  t  are  evaluated  in  the  same  manne  The  n  order  partial 
derivative  is  given  by  equation  (5): 

_  snF(x,8)-8n  ^(XjOj-s11  23f (x,o)  (5) 

3t  3t 


The  Laplace  transform  of  af (x,t) 

Sx 


„n-l 


5t 


n-i 


f (x,o) 


is  evaluated  as  follows: 


F(x,s)  = 


f (x, t)exp(-ts)dt 


dF<*;3)  * 
ax  ax 


f (x. t)exp(-ts)dt  = 


_a_ 

ax 


[f (x,t)exp(-ts)]dt 


or 


ax 


exP("ts)dt 


(6) 


A  comparison  between  equations  (1)  and  (6)  will  show  that 


rdffr.t)  j 

ax 


aF(x,S) 

ax 


(7) 
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Higher  order  derivatives  are  handled  in  an  analogous  manner. 

The  Laplace  transform  of  the  n^  order  partial  of  f(x,t)  with  x  is: 

=  (8) 

^  n  5  n 

9t  Sx 

Equations  (5)  and  (8)  can  be  used  to  obtain  the  Laplace  transform 
of  mixed  and  higher  order  partial  derivatives  as  shown  by  Haberman  (8)» 

Having  completed  the  frequency  analysis  of  the  partial 
differential  equation,  the" lumped  parameter  transfer  function  must 
tiu;n  be  inverted  into  the  time  domain  so  as  to  obtain  the  desired 
time  domain  representation.  It  was  mentioned  earlier  that  in  most 
cases  the  lumped  parameter  transfer  function  cannot  be  readily  inverted 
by  classical  methods.  Thus  the  problem  is  one  of  selecting. a  suitable 
functional  form  for  the  transfer  function  and  secondly  to  "fit"  the 
transfer  function  to  the  frequency  response  data  such  that  there  is 
good  agreement  in  the  time  domain  (2) . 

Little  can  be  said  about  selecting  a  suitable  transfer  function, 
other  than  it  must  be  readily  inverted  into  the  time  domain  and  should 
include  terms  characteristic  of  the  response  of  the  particular  system. 

To  obtain  the  best  statistical  agreement  in  the  time  domain, 


the  following  criteria  function  was  chosen  for  minimizations: 


t 

SSDRT  =  J8t  (g(x,t)  -h(x,t)}2  dt 
o 

where 


g(x,t) 
h(x, t) 


=  actual  time  function  of  the  system 
=  approximate  time  function  of  the  system 
=  settling  time 


(9) 
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Through  an  application  of  Parseval's  Theorem,  Schnelle  (3) 
has  sho^m  that  the  frequency  domain  equivalent  of  Equation  (9)  is: 

CO 

SSDRF  =  J  [G^  -  2GHCos(g<*0)-HI2]du>  (10) 

o 

where 

G  =  |f(x,u>)j  =  actual  magnitude  of  system  frequency  response 

H  =  ig(x,o))j  53  approximate  magnitude  of  system  frequency 

response 

<v  =  actual  phase  angle  of  system  frequency  response 

0  =  approximate  phase  angle  of  system  frequency  response 
Thus  to  obtain  the  best  time  domain  representation  by  the 
minimizations  of  Equation  (9),  the  integral  in  Equation  (10)  should  be 
minimized  when  working  with  data  in  the "frequency  domain. 


Development  of  Model  Equation 

A  relatively  simple  system  that  could  be  used  to  demonstrate  the 
application  of  this  method  to  distributive  parameter  systems  of  the 
convective  type  is  a  heat  exchanger  in  which  the  temperature  on  one 
side  of  the  tube  wall  is  constant. 

The  particular  system  chosen  consists  of  a  one  tube-pass,  one 
shell-pass  steam-hydrocarbon  heat  exchanger  with  the  condensing  steam 
on  the  shell  side.  (5) 

The  heat  duty  on  the  exchanger  will  be  sensitive" to  changes  in 
the  con,. .using  pressure  of  the  steam,  hydrocarbon  flow  rate  and 
hyc!  .-..mon  inlet  temperature. 

A  transfer  function  relating  the  liquid  hydrocarbon  outlet 
temperature  to  condensing  temperature  of  the  steam  will  be  developed 
for  the  particular  system  just  described  and  pictured  below. 


Figure  1.  Pictorial  representation  of  convective  system  under  study. 
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The  interactions  between  the  shell  wall  capacity  and  tube  wall 
capacity  are  neglected.  Also  the  lag  in  the  steam  pressure  to  the 
valve  opening  due  to  the  capacitance  of  the  shell  wall  and  vapors 
within  the  shell  is  not  considered.  In  order  to  maintain  a  constant 
tube  wall  temperature  throughout  the  length  of  the  exchanger,  it  is 
necessary  to  assume  that  the  tube  wall  has  no  resistance  and  negli¬ 
gible  capacitance  in  the  condensate  film  when  compared  with  the 
capacity  of  the  liquid  hydrocarbon  film.  The  more  common  assumptions 
of  constant  physical  properties,  no  backmixing  and  constant  liquid 
hydrocarbon  velocity  are  also  employed. 

By  making  use  of  these  assumptions,  the.  resulting  convective 
equation  is: 


Tl  dT(x, t)  Tlv  dT(x,t) 

dt  3x 


T  (t)  -  T(x, t) 
w 


(11) 


where 


A  pC. 

11  ■  riH  <1/sec) 

i  is 


(12) 


and  the  wall  equation: 


T? 


dT  (t) 


w 


where 


dt 


+  (T3+l)Tw(t)  =  Ta+T3T(x, t) 


(13) 


T2  = 


M  C 
w  w 

h  A 

n  SO 


(l/sec) 


(14) 


T3 


h  A  , 
1  si 

h  A 
o  so 


(ditnensionleae) 


(15) 


Equations  (11)  and  (13)  complete  the  mathematical  description  of 
the  system  under  study.  These  equations  are  coupled  through  the  wall 
temperature  indicating  that  a  simultaneous  solution  is  necessary  if 
the  system  is  to  be  solved. 

Frequency  Analysis 

It  is  apparent  that  the  simultaneous  solution  of  Equations  (11) 
and  (13)  via  classical  techniques  would  be  very  difficult  to  implement. 

An  alternative  method  for  handling  these  simultaneous  equations 
is  by  transformation  into  the  frequency  domain.  The  transformation 
from  the  time  domain  into  the  frequency  domain  reduces  the  wall  equation 
to  an  algebraic  equation  and  the  convective  equation  to  an  ordinary 
differential  equation. 

It  is  convenient  to  define  new  variables  representing  deviations 
from  the  steady  state  value  along  any  point  within  the  exchanger  before 
transformation  of  the  equations. 

6  (x,t)  «  T(x,t)-T  (x)  (16) 

88 

e  (t)  -  T  -(T  )ss  (17) 

w  w  w 

e  (t)  =  T  -(T  )ss  (18) 

8  So 

After  substituting  above  variables  in  Equations  (11)  and  (13) 
and  cancelling  the  steady  state  terms,  the  following  equations  result: 

vti  c  e  -e  (i9) 

dx  St  w 

de 

T2  ~  +  (T3+l)0w(s)  +  T30(x,s)  (20) 


9 


(28) 


Equation  (28)  is  a  transfer  function  relating  liquid  hydrocarbon 
temperature  at  any  point  within  the  exchanger  to  the  temperature  of  the 
condensing  steam.  To  obtain  the  desired  transfer  function,  Lc  is  only 
necessary  to  lump  the  distributive  liquid  hydrocarbon  temperature  at 
the  outlet  of  the  exchanger.  The  desired  transfer  function  relating 
the  lumped  liquid  hydrocarbon  temperature  to  changes  in  steam 
temperature  is: 


v«3 f-'-Mv#.,  ^5"- f  «S”S'-!.~  j-'-gfa-usKv ' 


&  r 


Computational  Procedure 


Upon  completing  the  development  of  the  model  equation  and  its 
corresponding  frequency  solution,  consideration  must  next  be  given  to 
the  computational  techniques  that  will  enable  their  nuaeeidal  solution. 
The  computational  effort  was, twofold; . one  aimed  at  the  solution  of 
the  model  equation  and  the  other  aimed  at  determining  an  acceptable 
transfer  function  to  approximate  the  frequency  solution  of  tehe  model 
equation. 

Comparison  of  the  time  domain  response  of  the  approximate 
transfer  functions  to  the  numerical  solution  of  the  partial  differential 
equations  then  provided  a  means  of  evaluating  the  postulated  transfer 
function. 

By  a  linear  transformation  of  the  existing  variables  and 
coordinates,  equations  (11)  and  (13)  can  be  expressed  in  a  more  general 
form  that  at  the  same  time  is  conducive  to  numerical  calculations. 

Defining  a  new  set  of  temperatures  of  magnitude  dictated  by  the 


arbitrary  parameter  Tq. 


U(y,  Y)  =  ^-t>-T(0»t> 


(30) 


a 


TlJ  \  J  —  ■> 

V(y,Y)  -  ~ — j - 

(31) 

a 

T  (t)-T(o,t) 

W(y,Y)  =■  - 

(32) 

Defining  now  the  new  special  and  time  coordinates  through  the 


usage  of  arbitrary  parameters  x  and  t  . 

fl  & 
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With  the  new  set  of  temperatures  a:  I  coordinates,  the  convective 
and  wall  equations  become: 


n  +  Tlv  .  V(»,u(v  „v 

T  BY  x  By  V'‘Y'  u'y»*' 

a  a 


?^d?^  +  fls+Wvco-Wy.t)  ■  >t(y) 

a 


By  choosing  T  to  represent  the  algebraic  difference  between  the 
a 

initial  steam  temperature  and  the  inlet  liquid  hydrocarbon  temperature, 

the  temperatures  under  consideration  can  be  made  dimensionless.  To 

normalize  the  spacial  coordinates,  x  is  chosen  to  be  the  length  of  the 

a 

exchanger.  By  assigning  to  t  the  numerical  value  of  the  time  dependent 

a 

T2,  the  time  coordinate  is  made  dimensionless  and  at  the  same  time  a 
slight  simplification  of  the  wall  equation  is  accomplished.  Equations 
(36),  (37)  and  (38)  specify  the  value  of  the  parameters  just  discussed. 


T  (o)-T(o,t) 

D 


t  -  T2  ( 

a 

By  using  the  above  values  assigned  to  T^,  xfl  and  tfl,  the  con¬ 


vective  and  wall  equations  become: 


li  MbXL+  ELIESliXL  =  v(Y)“U(y,v') 

T2  ay  L  By  ; 


+  (T4+l)V(Y)-l3U(y,Y)  -  W(Y) 
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i»,'5v«!?»'M55 


%  ,-^ss  >  “W’y*;  >t<^-*  ?/^nv^73 


An  abbreviated  from  of  the  above  equation  results  by  defining. 


C  a  — *• 

v/.  m/) 


T3+1 

(39) 

-T3 

(40) 

Tl 

T2 

(41) 

Tlv 

L 

(42) 

This  form  being 


C4  '~U~^^'  +  B  -  V(T)-U(y,t) 


+  C0V(Y)  *  C,U(y,Y>  -  W(Y) 


(43) 


(44) 


2  3 

A  numerical  solution  of  Equations  (43)  and  (44)  was  obtained 
using  central  difference  analogs  for  the  partial  derivatives. 

The  steady  state  solution  of  the  model  equation  was  also  obtained 
directly  from  equations  (43)  and  (44)  by  eliminating  the  terms  containing 
the  deeivative  of  the  dimensionless  time  coordinate.  The  resulting 
equations  possess  an  analytical  solution  given  below: 

U(y,o)  «  W(o)(l-exp(-G1y))  (45) 

where 

C, 


1 

BC, 


W(o)  “  C  V(x)  +  C3U(y,o)  (40 

This  initial  distribution  is  required  for  the  solution  of  the 
partial  differential  equations. 


m  & 
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»y  yw.  Iff.  > 


e  (s) 
s 


Steam 

Temperature 


|  $(**«) 

Gain  exp(-t  $) 

- -  Liquid  hydrocarbon 

(t^s+1)  temperature  at  outlet 


a)  First  order  lag  with  dead  timeJ..G^(s) 


Gain  exp(-tns)  w. 

Steam  '  “  “  Liquid  hydrocarbon 

Temperature  (tjS+l) (t2s+1)  temperature  at  outlet 


b)  Second  order  lag  with  dead  time...G2(s) 


(t^s+1) (t2s+1) 


+  0  (L ,  s ) 


Temperature 


Sgain  exp(-t  s) 
_ o 

(t^s+1) (t2s+1) 


~  Liquid  hydro-? 
carbon  tempera¬ 
ture  at  outlet 


Figure  2. 


c)  Combination  of  a  second  order  lag  and  a  second 
order  lag  with  dead  time...G^(s) 

Postulated  transfer  functions  used  in  analytical  lumped 
parameter  approximation  technique 


As  mentioned  earlier  in  this  section,  the  other  computational 
effort  was  aimed  at  obtaining  an  acceptable  transfer  function  that 
approximates  the  frequency  solution  of  the  model  equation  (29).  By 
using  equation  (29)  to  generate  frequency  data,  non-linear  regression 
analysis  was  carried  out  on  the  postulated  transfer  functions  and 
the  optimum  parameters  were  thus  obtained. 

A  program  based  on  equation  (29)  was  written  to  bbtain  the  - 
necessary  frequency  data  for  the  regression  analysis.  A  block 
diagram  representation  of  each  of  the  postulated  transfer  functions 
is  shown  in  Figure  2. 

For  comparative  purposes,  the  regression  analysis  was  carried 
out  on  two  objective  functions.  The  first  utilized  the  standard 
practice  of  obtaining  the  best  parameter  fit  in  the  frequency 
domain  by  minimizing 

CO 

PI  »  |  {(G  -  H)2  +  (9  -  0) 2}dci)  (47) 

This  will  be  referred  to  as  frequency  domain  minimization. 

The  second  objective  function  involved  minimizing 

|  {g(x,t)  -  h(x,t)}2  dt 
o 

using  frequency  domain  data  by  minimizing 

PI  -  J  {G2  -  2GH  Cos  (8  -  0)  +  H2}  du)  (48) 

o 

as  previously  discussed.  This  second  PI  will  be  referred  to  as  time 


domain  minimization. 


The  rotational  discrimination  program  developed  by  Farias  and 
Law  (i*)  was  used  to  carry  out  the  non-linear  regressions.  The  program 
is  based  In  a  modified  accelerated  gradient, technique,  and  is 
written  in  a  self-calling  subroutine  structure. 

For  each  transfer  function,  eoveral  regressions  were  carried 
out  in  both  the  frequency  and  time  domain.  Each  regression  was  started 
at  a  different  initial  point  with  tlie  hope  of  obtaining  the  best  set 
of  parameters. 

After  obtaining  the  beat  possible  set  of  parameters  for  each  of 
the  postulated  r.ransfer  functions,  their  frequency  response  was 
compered  with  that  of  equation  (29).  This  comparison  allowed  the 
determination  of  the  transfer  function  that  best  represented  the 
frequency  solution  of  the  model  equation.  Upon  selecting  the  postulated 
transfer  function  that  best  approximated  equation  (29),  the  final 
check  to  the  analytical  lumped  parameter  approximation  technique 
was  accomplished  by  inverting  this  transfer  function  into  the  time 
domain  and  comparing  it  with  the  numerical  solution  of  the  model 
equation. 

To  compare  time-domain  regressions  against  frequency  regressions, 
the  postulated  transfer  function  was  inverted  into  the  time  domain 
using  the  optimum  set  of  parameters  from  both  the  frequency  and  time- 
domain  regressions.  For  each  transfer  function  its  time  domain  repre¬ 
sentation  is  compared  against  the  numerics!  solution  of  the  model 
equation  by  first  using  the  frequency  regressed  parameters  and  then 
the  time-domain  regressed  parameters. 


Result  8 


Numerical  Solution  of  Model  Equation;  The  system  was  excited 
with  a  10° F  step  change  in  the  steam  temperature.  The  numerical 
integration  was  performed  with  several  AY/ AT  ratios,  with  a  ratio 
of  0.05  (0.005/0,10)  giving  the  best  results. 

Although  no  stability  problems  were  encountered,  a  slight  dis¬ 
crepancy  (in  the  fifth  significant  figure)  was  observed  between  the 
steady  state  values  of  the  liquid  hydrocarbon  temperature  calculated 
from  the  numerical  solution  and  those  calculated  analytically  from 
the  steady  state  equations.  This  discrepancy  was  due  to  the  capaci¬ 
tance  of  the  hydrocarbon  film  causing  the  tube  wall  temperature  to 
show  slight  distributive  tendencies.  Table  I  shows  the  steady  states 
calculated  from  the  numerical  solution  and  those  calculated  from  the 
steady  state  equations. 

The  time  domain  response  of  the  liquid  hydrocarbon  temperature 
at  the  outlet  of  the  exchanger  is  shown  in  Figure  3. 

Frequency  Analysis  of  Model  Equation:  The  frequency  response  of 
tne  model  equation  is  shown  in  Figure  4.  Note  the  resonance  peaks 
characteristic  of  this  system  occurred  in  the  medium  and  high 
frequencies  ((bs0.55  rads/sec).  From  this  data,  thirty-five  points 
were  selected  for  the  non-linear  regressions. 

Based  on  the  data  selected  from  the  frequency  response  of  the 
model  equation,  several  starting  points  were  used  for  each  of  the 
postulated  transfer  functions  with  the  hope  of  obtaining  a  good  estimate 
of  the  global  optimum.  Tables  II  and  III  show  the  unimodel  character¬ 
istics  of  the  response  surface  resulting  from  the  frequency-domain 
regressions  of  G^(s)  and  62(8).  These  two  tables  show  the  local 
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TABLE  I 


Actual 

STEADY  STATES  OF  THE  DISTRIBUTED 

LIQUID  HYDROCARBON  TEMPERATURE 

Initial  Steady  State  Final  Sto 

ady  State 

Position 

Numerical 

Steady  State 

Numerical 

Steady  State 

(Ft) 

Solution 

Equation 

Solution 

Ecmation 

(F) 

(F) 

<“F) 

CF) 

0.0 

70.000 

70.000 

70.000 

70.000 

0.4 

74.444 

74.434 

74.596 

74.589 

O.S 

78.818 

78.800 

79.119 

79.107 

1.2 

83.123 

33.098 

83.573 

83.556 

1.6 

87.361 

87.331' 

87.956 

87.934 

2.0 

91.532 

91.497 

92.271 

92.246 

2.4 

95.638 

95.599 

96.519 

96.491 

2.8 

99.679 

99.638 

10f.700 

100.671 

3.2 

103.656 

103.615 

104.815 

104.786 

3.6 

107.571 

107.530 

108.867 

108.837 

4.0 

111.425 

111.384 

112.855 

112.826 

4.4 

115.218 

115.179 

115.781 

116.753 

4.8 

118.951 

118.915 

120.645 

120.620 

5.2 

122.626 

122.594 

124.449 

124.427 

5.6 

126.243 

126.216 

128.193 

128.175 

6.0 

129.803 

129.782 

135.879 

131.865 

6.4 

133.307 

133.292 

135,508 

135.498 

6.8 

136.757 

136.749 

139.379 

139.074 

7.2 

140.152 

140.152 

142  .395 

142.596 

7.6 

143.493 

143.502 

146.056 

146.063 

8.0 

146.782 

146.801 

149,463' 

149.477 

i, 


1.00  0.0822956  _2 

10.00  3:28888  2.405361x10  Time 

10.00  3.28684 


Table  III  (Cont'd) 
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optima  that  were  encountered  during  the  time-domain  regressions,  and 
hence  a  considerable  disadvantage  over  the  frequency-domain  regressions 
of  G^(s)  and  GgCs).  A  considerable  number  of  local  optima  were 
encountered  during  the  frequency  and  time-domain  regressions  on 
G^Cs)  as  shown  in  Table  IV.  It  appears  that  the  simultaneous 
optimization  of  the  dead  time  and  time  lags  tend  to  introduce  these 
local  optima  in  the  response  surface.  The  simultaneous  regressions 
are  inherent  to  time  domain  optimizations  and  also  were  carried  out 
in  the  frequency-domain  regressions  of  Gg(s).  Table  V  shows  the 
optimum  parameters  that  were  determined  for  both  frequency  and  time- 
domain  regressions  for  each  of  the  postulated  transfer  functions. 

The  frequency  response  of  G^(s)  for  both  frequency  and  time- 
domain  regressed  parameters  are  shown  in  Figure  5^  A  slight 

difference  in  the  responses  of  the  frequency  and  time-domain 
regressed  parameters  was  observed.  Both  responses  gave  poor 
representation  of  tht  frequency  response  of  the  model  equation 
through  the  high,  middle  and  part  of  the  low  frequencies. 

The  frequency  response  of  G2(s)  is  shown  in  Figures  6  and  7. 

No  substantial  difference  in  the  responses  of  the  frequency  and 
time-domain  regressed  parameters  was  observed.  The  frequency 
response  of  G2(s)  gave  a  better  frequency  approximation  in  the  low 
and  part  of  the  middle  frequencies,  while  the  high  and  part  of  the 
middle  frequencies  still  remain  poorly  represented. 

The  frequency  response  of  Gg(s)  is  shown  in  Figure  8.  In  this 
case  too,  no  appreciable  difference  in  the  response  of  the  frequency 
and  time-domain  regressed  parameters  was  observed.  G^Cs)  gave  an 
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11.00  11.4299 

37.00  36.6997  1.483555x10  Freq 

0.33  0.334694 


Tigure  5a.  Frequency  response  of  G^(s)  obtained  in  the  analytical 
lumped  parameter  approximation  technique. 
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Figure  8b,  Frequency  response  of  C^(s)  obtained  in  the  analytics 
lumped  parameter  apvn- .motion  technique. 
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ponse  of  G  (s)  obtained  in  the  analytical  lumped  parameter 


almost  perfect  frequency  representation  of  the  model  equation 
throughout  the  frequency  range  studied. 

Time  Domain  Analysis  of  Model  Equation:  The  time  domain 
representation  of  G, (s)  for  both  frequency  and  time-domain  regressed 
parameters  are  shown  in  Figure  9.  With  either  frequency  or  time-domain 
regressed  parameters,  this  transfer  function  gives  a  rather  poor 
representation  of  the  liquid  hydrocarbon  temperature  at  the  outlet  of 
the  exchanger  at  any  given  instant  of  time.  A  small  difference  in 
responses  between  the  frequency  and  time-domain  regressed  parameters 
giving  a  slightly  better  approximation  to  the  dynamics  of  the  liquid 
hydrocarbon  temperature . 

The  time  domain  response  of  62(8)  is  shown  in  Figure  10.  No 
difference  was  observed  between  the  responses  of  the  frequency  and 
time-domain  regressed  parameters  in  this  case.  A  better  approximation 
of  the  liquid  hydrocarbon  temperature  dynamics  was  achieved  with  this 
transfer  function  particularly  during  the  first  part  of  the  solution 
time. 

The  time  domain  representation  of  G^s)  is  shown  in  Figure  11. 

In  this  case  too,  no  difference  was  observed  between  the  responses  of 
the  frequency  and  time-domain  regressed  parameters.  This  transfer 
function  provided  an  excellent  description  of  the  outlet  liquid 
hydrocarbon  temperature  dynamics  throughout  the  solution  time. 

Conclusions 

Excellent  representation  of  the  time  response  of  the  liquid 
hydrocarbon  temperature  at  the  outlet  of  the  exchanger  was  obtained 
with  the  combination  of  second  order  lag  and  second  order  lag  with 


dead  time  transfer  function.  Another  advantage  associated  with  this 
transfer  function  is  that  it  only  contains  three  parameters  and  can 
be  easily  implemented  in  analog  simulation  studies. 

The  lack  of  improvement  over  frequency  regressions  coupled  to 
the  additional  computational  effort,  makes  time-domain  regressions 
impractical' for  convective  systems. 
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NOMENCLATURE 


2 

A.  flow  area  of  exchanger  tubes,  ft 

ic 

2 

inside  surface  area  of  exchanger  tubes,  ft  /ft 

2 

Ago  outside  surface  area  of  exchanger  tubes,  ft ‘/ft 

Heat  capacity  of  liquid  hydrocarbon  Btu/lb  °F 
f(x,t)  8 pace  dependent  variable 

F(x,s)  frequency  equivalent  of  f(x,t) 

G^(s)  first  order  lag  with  dead  time  transfer  function 
G2(s)  second  order  lag  with  dead  time  transfer  function 

G^(s)  combination  of  second  order  lag  and  second  order  lag  with 

dead  time  transfer  function 
Gain  steady  state  gain  of  system 

2 

h^,  inside  heat  transfer . coefficient,  Btu/hr(°F)ft 

2 

hQ  outside  heat  transfer  coefficient,  Btu/hr(°F)ft 

L  exchanger  length,  ft 

Mw  linear  density  of  exchanger  tubes,  Ib/ft 

P  8 team  pressure,  psig 

Sgain  1-Gain 

t  time,  seconds 

t^  duration  of  input  pulse,  seconds 

t  dead  time,  seconds 

o 

t  duration  of  output  pulse,  seconds 

T(x,t)  distributed  liquid  hydrocarbon  temperature,  °F 

T  (t)  8 team  temperature,  °F 

s 

Tw(t)  tube  wall  temperature,  °F 

T  (x)  steady  state  value  of  the  distributed  liquid  hydrocarbon 
88 

temperature,  °F 


(T  )  • 

W  88 

(T  ) 

8  88 

U(y,Y) 


v 

V(Y) 

W(Y) 


x 


y 


steady  state  value  of  the  tube  wall  temperature,  °F 

steady  state  value  of  the  steam  temperature,  °F 

distributed  liquid  hydrocarbon  temperature,  dimensionless 

liquid  hydrocarbon  velocity  within  exchanger  tabes,  ft/sec 

tube  wall  temperature,  dimensionless 

steam  temperature,  dimensionless 

distance  along  exchanger,  ft 

distance  along  exchanger,  dimensionless 

unit  step  applied  at  time  »  t 

o 


Greek 

Y 

U) 

P 

0(x,t) 


88<t> 

e„<t> 


Q(x,s) 

«.<■) 

v> 


time,  dimensionless 
frequency,  radians /second 

3 

liquid  hydrocarbon  density,  lb/ ft 

distributed  liquid  hydrocarbon  temperature  representing 
deviations  from  steady  state  value,  °F 
steam  temperature  representing  deviations  from  steady 
state  value,  °F 

wall  temperature  representing  deviations  from  steady 

state  value,  °F 

frequency  equivalent  of  0(x,t) 

frequency  equivalent  of  9  (t) 

8 

frequency  equivalent  of  @w(t) 
time  constant,  seconds 
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